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Abstract 

We solve the Gauss law and the corresponding Mandelstam constraints in the loop Hilbert space 
Ti.^ using the prepotential formulation of (d+ 1) dimensional SU(2) lattice gauge theory. The resulting 
orthonormal and complete loop basis, explicitly constructed in terms of the d{2d — 1) prepotential 
■ intertwining operators, is used to transcribe the gauge dynamics directly in Ti^ without any redundant 

' gauge and loop degrees of freedom. Using generalized Wigner-Eckart theorem and Biedenharn -Elliot 

identity in Ti^ , we show that the above loop dynamics for pure SU(2) lattice gauge theory in arbitrary 



o 

(N 



dimension, is given by real and symmetric 3nj coefficients of the second kind (e.g., n=6, 10 for d=2, 3 
^ ■ respectively). The corresponding "ribbon diagrams" representing SU(2) loop dynamics are constructed. 

. The prepotential techniques are trivially extended to include fundamental matter fields leading to a 

. description in terms of loops and strings. The SU(N) gauge group is briefly discussed. 

m ■ 

> : 

■ 1 Introduction 
O ■ 

o ; 

, The idea that gauge theories should be formulated completely in terms of loops in space carrying electric 

, fluxes is quite old, appealing and has long history [TJ [51 [31 [H [51 [B] . In the context of electrodynamics, 

, Yang [S] has emphasized the importance of path dependent "non-integrable phase factors" carrying 

, electric fluxes to provide a complete description of all the quantum effects. In the context of quantum 

I ' chromodynamics the loop formulation, without any colored gluon or colored quark degrees of freedom, 

\ is further expected to provide a better framework to analyze non-perturbative low energy issues like 

'"^ ' color confinement. In the context of gravity, the relatively recent Hamiltonian formulation of quantum 

^ ' gravity in terms of SU(2) connections has also been reformulated in terms of loops leading to loop 

' quantum gravity [7j. Therefore, the loop formulation of gauge theories may eventually provide a 

' common geometrical platform to understand all interactions in nature. In this work, we analyze some 
of the basic kinematical as well as dynamical issues involved in the Hamiltonian formulation of lattice 
gauge theories in terms of loops. 

The standard construction of the loop states for pure SU(N) gauge theory on lattice [S] consists of 
considering the set of all oriented loops T and constructing the corresponding gauge invariant Wilson 
loop operators = TrU^ for every 7 € F. Acting on the strong coupling vacuum, all possible gauge 
invariant operators of the form W-y^W^^....Wy^ create all possible gauge invariant states associated 
with the corresponding loops 71 72... 7m. These loop states are manifestly gauge invariant, geometrical 
and form a basis, usually known as Wilson loop basis (see section 2.1). However, a serious problem with 
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this Wilson loop basis is that it again over-describes gauge theory. This time, the over description 
is because loops of all shapes and sizes have to be included in constructing the above loop basis. 
Therefore, one is again confronted with too many redundant (but now SU(N) gauge invariant) loop 
degrees of freedom (see section 2.1). The Mandelstam constraints Ij amongst the various loop states 
express this over-completeness of the Wilson loop basis (see section 2.1). As these constraints represent 
the linear dependence of the gauge invariant states associated with loops of all sizes, they are highly 
non-local and hence difficult to solve (see section 2.1). In the strong coupling {g — > oo) limit, the loops 
are small and they carry small electric fluxes [5]. Therefore, the Mandelstam constraints can be easily 
solved by using Gram-Schmidt orthogonalization procedure amongst the small number of loop states 
considered within this (g — > oo) limit. However, in the continuum (g 0) limit, as opposed to the 
strong coupling limit, large loops carrying large electric fluxes will be important |15j . Therefore, the 
problem of over-completeness of the Wilson loop basis will become more and more acute as we remove 
the lattice cut-off and approach the continuum limit. This over-completeness, in turn, will result in 
rapid proliferation of spurious zero modes of the Hamiltonian in the Wilson loop basis. Therefore, the 
initial problem in loop formulation is to solve the Mandelstam constraints exactly before addressing 
any dynamical issue. Infact, as stated by Gambini and Pullin in [sj^: "The proliferation of loops when 
one considers larger lattices and higher dimensions completely washes out the advantages provided by 
the (loop) formalism." . 

The motivation and purpose of the present work is to systematically develop ideas and techniques 
to reformulate lattice gauge theories in loop space without any spurious loop degrees of freedom. 
We solve SU(2) Mandelstam constraints leading to an orthonornial loop basis which is complete and 
characterized exactly by 3{d — 1) angular momentum quantum numbers per lattice site [lOj . Further, 
we show that the most economical description of the pure (d+l) SU(2) lattice gauge theory dynamics, 
involving only the relevant loop degrees of freedom, is given by 3nj coefficients of the second kind and 
therefore it is highly geometrical. 

We will work within the prepotential formulation OHU] of SU(2) lattice gauge theory Hamiltonian [S]. 
The prepotential approach has extended SU{2) (g) U{1) gauge invariance. The prepotential operators 
are SU(2) harmonic oscillator doublets attached to the initial and the final points of every link. They 
are connected by the U(l) gauge invariance mentioned above (see section 2.2). The advantage of 
prepotential approach is that the non-local and over-complete Wilson loop basis can be constructed 
and studied locally in terms of the SU(2) invariant d{2d — 1) "prepotential intertwining operators" at 
every lattice site. This local description of the Wilson loop basis is characterized by 2d{d — 1) (integer) 
intertwining quantum number per lattice site. Further, the Mandelstam constraints, which appear 
highly non-local in terms of the link operators, become local in terms of the prepotential operators 
(see section 2.3.1). This enables us to solve them explicitly using simple group theoretical ideas 
[TT1[T2]. The novel U(l) gauge invariance of the prepotential formulation then enables us to explicitly 
construct an orthonormal and complete loop basis in the entire loop Hilbert space Ti.^ on the lattice. 
As expected, this orthonormal basis is characterized by 3{d — 1) (as opposed to 2d{d — 1) quantum 
numbers for the Wilson loop basis) gauge invariant quantum numbers per lattice site. Having solved 
all the constraints, we use the generalized Wigner-Eckart theorem and Biedenharn-EUiot identities in 
2 Chapter 12, page 303-304 
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the resulting orthonormal loop basis to compute the loop dynamics. Our approach also enables us to 
compute loop dynamics locally site by site. The final results at different lattice site are then glued 
together through the U(l) gauge invariance. We show that in d+1 dimension the SU(2) loop dynamics 

in the (/, J) plane where /, J = 1 d and I < J plane is given by 3nj coefficients of the second kind 

where n = 2 [2 + d + ( J — /) — 6i,i — Sj,d]- The matrix elements of the Hamiltonian in the above loop 
basis are found to be real and symmetric. The 3nj coefficients of second kind, and therefore the loop 
dynamics, are graphically represented by the ribbon diagrams. 

The plan of the paper is as follows. In the first half (section 2) we discuss the kinematical issues 
and in the second half (section 3) we discuss the dynamical issues. In both these sections, the explicit 
computations are done in d = 2. This keeps the discussions simple and also illustrates all the essential 
ideas and techniques involving prepotentials. Their generalization to arbitrary d dimension is obvious 
and done next. As the complications caused by over-completeness of the loop basis or equivalently 
the Mandelstam constraints have been major obstacles in the loop approach to gauge theories, we 
first review them on lattice in section 2.1. In section 2.2, we briefly discuss the SU{2) (g) U{1) gauge 
invariant prepotential approach |^ which enables us to cast the Mandelstam constraints in a simple 
local form. In section 2.3, we solve the Mandelstam constraints and give all possible orthonormal loop 
state solutions pTOl in terms of the d{2d — 1) prepotential intertwining operators. In section 2.4, we 
discuss inclusion of matter fields leading to a gauge theory description in terms of loops and strings. 
In section 3.1 and 3.2, we compute the matrix elements of the Hamiltonian in the above loop basis and 
discuss the ribbon diagrams representing these amplitudes. In section 4, SU(N) gauge group is briefly 
discussed. The techniques used in constructing the orthonormal loop state basis are given in appendix 
A. The technical details involved in computing loop dynamics are given in appendix B. 



2 The kinematical issues 



This section is devoted to the kinematical issues involving the loop states in pure SU(2) lattice gauge 
theory. We first discuss the Mandelstam constraints in terms of the original lattice link operators. They 
look highly non-local. We then cast them in their local form at every lattice site n in terms of the SU(2) 
invariant prepotential intertwining operators at n. Next, we convert the problem of solving these local 
Mandelstam constraints to the problem of finding common (orthonormal) eigenvectors of a complete 
set of commuting observables (C.S.C.O) containing (4d — 3) angular momentum operators at n. In 
particular, we show that the states related by the Mandelstam constraints at n are degenerate with 
respect to a subset of operators belonging to the above C.S.C.O.. Therefore, the common eigenvectors 
of the C.S.C.O at n, which are specific linear combinations of the states related by the Mandelstam 
constraints, lead us to an orthonormal local basis at n. The final orthonormal loop basis over the entire 
lattice, characterized by 3(d — 1) angular momentum quantum numbers per lattice site, is obtained 
by weaving or gluing these local orthonormal basis at different lattice sites according to the additional 
U(l) Gauss law associated with the prepotential formulation. 
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2.1 Mandelstam constraints on lattice 

On lattice the number of gauge invariant degrees of freedom (Af) is given by the dimension of the 
quotient space ®iinksSU{2)/ ®sites SU{2). Thus for a d-dimcnsional periodic lattice with n'^ sites and 
dn'' links: 

Af ^3{d-l)n''. (1) 

Therefore, a complete description of the SU(2) gauge invariant physical Hilbert space on a periodic 
lattice should require J\f quantum numbers or equivalently 3{d — 1) quantum numbers per lattice 
site. However, the Wilson loop basis is characterized by = 2d{d — linking or intertwining 
quantum numbers (see section 2.3). Thus there are M. — N = (2d^ — Sd + 3) n'^ redundant degrees of 
freedom. This clearly establishes that the Wilson loop basis is over-complete. Further, the redundant 
degrees of freedom or degree of over-completeness increases as d^ for large d. Let us illustrate this 
over-completeness by giving the simplest example in d = 2. We consider two plaquettes A and B 
touching each other at a common lattice site n as shown in Figure ([T^). The corresponding Wilson 
loop operators satisfy: 

(TtUa) (TtUb) = Tr (UaUb) + Tr {UAUg') . (2) 

The relation Q is a trivial identity involving any two SU(2) matrices Ua and Ub- It can be checked 
by writing Ux = XqI + 'iJ2^=i^a.<7'^ where ct° are the Pauli matrices, Xo,Xa are real and satisfy 
Xq + Xf + X2 + = 1 . The generalisations of the identity ^ to the other classical Lie groups have 
been discussed in ^13j. We define the following three loop states: 

|7i) = (Tr;7A)(TrC/s)|0), I72) = Tr (t/^f/^i) |0), I73) = Tr (C/^[/b) |0), (3) 

The identity ^ implies the simplest Mandelstam constraints: 

I71) = |72> + |73)- (4) 



(a) 



(b) 
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n 

B 

> 


— :> 



(c) 



Figure 1: The graphical representation of the simplest Mandelstam constraint ^ with Na — Nb — 1- 
The thick lines are for the later comparison of the same in the prepotential formulation. 
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Thus we see that the three loop states I71), I72) and I73) are hnearly dependent. To appreciate the 
problem further, let us consider most general loop states involving only these two plaquettes A and B: 



\Na,Nb) = (TrUAf-'iTTUBf^lO) 

= (TrC/A)'^^"'(TrC/B)^«-i(TrC/AC/B + TrUAUg^)\0) 
= {TTUAf-'-^{TTUB)^^-\TTUAUB + TrUAUg^flO) 



= (Tr[/A)^-*-^—(Tr[/B)^«-^—(TrC/At/B + Tr[/A[/3i)^"""|0) (5) 

where Na, Nb are two arbitrary integers representing the angular momentum fluxes over A and B 
and Afinin=Minimum(Af^, A^b). Thus, given the loop state \Na,Nb), we have produced {2Nnun + 1) 
distinct but linearly dependent Wilson loop states contained in the relations ([5]). Therefore, even in 
d = 2, for the simplest loop states over the two plaquettes, the number of linearly dependent loop 
states increases with the SU(2) flux value A'min. Note that we have to include states with Nmin taking 
arbitrarily large values. Next one can imagine extending the picture in Figure ([T]) by adding more 
plaquettes leading to further constraintifl. In general, one has to address the problem of finding a 
complete set of linearly independent loop states amongst the loop states of all shapes, sizes carrying 
arbitrary fluxes and touching/crossing one another at arbitrary number of lattice sites. Further, it is 
clear that this problem becomes more and more difficult as the dimension d increases. This is the 
reason why the Mandelstam constraints have been notorious and major obstacle in the development 
of the loop formulation of gauge theories . The solutions of the Mandelstam constraints have been 
discussed at various levels of approximations and difficulties in the past. In [H [H] an approximate 
loop cluster method in 2+1 dimensions is developed and the Schrodinger equation is written down as 
difference equations in these cluster coordinates. However, approximations involved are uncontrolled 
[5] to enable us to go to 3+1 dimensions. Motivated by strong coupling expansions, various methods 
involving cluster of loops and their truncations to define Hamiltonian eigenvalue problem have been 
developed to go beyond strong coupling region [5D]. In ^Ij the Mandelstam constraints are solved 
and eigenvalues equations are analyzed on computer using small lattices and small loops. In [231 112 
the Mandelstam constraints are solved classically on a finite periodic d=2, 3 lattice leading to an 
independent loop configuration space. However, the issues like quantization of these loop variables 
and setting up the corresponding Schrodinger equation are not clear [53]. We now briefly review the 
various operators in the Kogut-Susskind formulation and define the prepotential operators [9 which 
enable us to solve the Mandelstam constraints exactly and explicitly in arbitrary dimensions. 

2.2 The prepotential operators 

The kinematical variables involved in Kogut and Susskind Hamiltonian formulation [5J of lattice gauge 
theories describe SU(2) rigid rotators attached to every link (n, i) of the lattice. The kinematical 

^In |6] it is shown that the Mandelstam constraints constitute sufficient algebraic conditions on Wilson loop variables 
to allow reconstruction of the corresponding gauge potentials. 
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, \ ( Uii{n,i) Ui2{n,i) \ 
variables are: a) SU(2) link operators U{n,i) = \ , describing the orientation 

\ U2i{n,i) U22{n,i) J 

of the body fixed frame of the rigid body from the space fixed frame, b) the electric fields E1{n, i) and 
E'^in + which are the components of the angular momentum in the body fixed and space fixed 
frames respectively. This description is shown in Figure ^^). The link operators U{n,i) satisfy the 
SU(2) conditions: 

U{n,i)U''{n,i) = U''{n,i)U{n,i)=I, \U{n,i)\ = l. (6) 
Above X is 2 X 2 identity matrix and \U\ = dctU. The rigid body commutation relations arJ3 [S]: 

[El{n,i),U{n,i)]^-(^^^ U{n,t) ^ [El{n,z), El{n,z)] = le^'^Elin^i), 

[E%in + z,i),U{n,t)]^U{n,z)(^^^ => £;^(n, »)] = ze^'^E^j.in^t) (7) 

In ([71), cr^ (a=l,2,3) are the Pauli matrices and the operators on different hnks commute. The angular 



rpa 



momentum algebras on the r.h.s. of ([7|) follows from the Jacobi identities: 

cyclic permutations = 0. The link operators commute amongst themselves: [Uap{n,i),U-ys{m, j) 

Uaf3{n,i),U}^g{'m, j) = 0. Further, E1{n,i) and E'^{n + i,i), being the body, space fixed compo- 
nents of the angular momentum operator of the rigid rotator on the link (n,i), mutually commute 
[El{n,i),E\{n + =0. For the same reason, they satisfy the kinematical constraint: 



3 3 

E 

a— 1 a— 1 



i?£(n,*)ii;£(n,z)=^i?^,(n + z,z)i?^(n + *,z) = i?2(n,*), V(n, z) (8) 



ensuring that their magnitudes are equal. The SU(2) gauge transformations correspond to separately 
rotating the body, space fixed frames of the rigid rotator [5]: 

EL(R){n,i)^ K{n)EL(R){n,i)h.\n), U{n,i) ^ K{n)U{n,i)K\n + i). (9) 

In E]^(^fi) = X)a=i ^L{R)^'^ ^'^'^ A(n) is rotation matrix in the fundamental representation of SU(2) 
at lattice site n. The commutation relations ^ along with the gauge transformations ([9]) imply that 
the generators of SU(2) gauge transformations at any lattice site n are: 

3 

C'^{n)^Y.{^LM+El{n,i)), yn,a. (10) 

1=1 

The corresponding Gauss law constraints are C(n) — 0. The commutation relations ([7]) and the con- 
straints ((8|) imply that a complete set of commuting observables on every link (n, j) are: E^ (n^i), Ef^ [n^i) 
E^{n + where E^i^j^s^{n,i) = E'^^^in^i) are the third components of the angular momenta. The 
corresponding orthonormal basis is denoted by |j(n, i), TO(n, i), TO(n -I- i, i)) where j(n,i), m{n,i), 
and m(n -|- are the eigenvalues of the above 3 mutually commuting operators respectively. We 
now define the prepotential operators through the Jordan-Schwinger representation of the angular 
momentum algebra [17j : 

El{n,i) = iat(n,i)cr°a(n,i), E%{n,i) = ^&^(n, f)cr°6(n, f). (11) 



*In [9], we had used E°'{n, i) = —E'^{n, i) and e"(ra, i) = E'^{n+i, i) resulting in an extra -ve signs in the commutation 
relations involving i?£ in 
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a^'.a = b«b 



E^n,i) 



U(n,i) 



E (n+i,i) 
R 



aV,i) 
a 



+ 

b (n+M) 



Figure 2: (a) The original {E1{n,i), E1^{n + i,i),Uap{n,i)) operators, (b) the new 
(00(71., j), 6ct(n + «, i)) operators on the Hnk {n,i). The electric fields EljEr transform in the 
adjoint representation. The prepotentials aa(ji,i),ba(n + i,i) transform as SU(2) fundamental 
doublets at lattice site n and n+i respectively. Therefore, we represent them by Young tableau □ at 
n and n+i respectively. 

The mapping ([TT|) corresponds to associating two doublets of harmonic oscillator prepotentials a]j(n, i) 
and 6jj(n + j, i) and their conjugates to the initial and the end points of the link {n, i) respectively. This 
assignment is shown in Figure ([Ih)- The canonical electric field or angular momentum commutation 
relations in ^ are satisfied provided the prepotentials on every link (n, i) satisfy the standard harmonic 
oscillator algebra: 



; ^p 



Saf3, 



5a0, [aq, a/3] = 0, [6q, 6/3] = 0. 



(12) 



The body and space fixed components of the angular momentum or electric fields mutually commute 
implying: [ac(,6jj] = [0^,6^] = 0. Note that the prepotential vacuum state |0) on the link (n,i) 
satisfying: aa{n,i)\0) = 0, ba{n,i)\Q) = is the strong coupling vacuum defined as: -E£(n,i)|0) = 
0, E'^(n,i)\Q) = 0. In [24] anti-commuting oscillators, instead of (fT2|) . are used in (fTTjl to treat QCD 
as quantum link models. We define the total number operators: 

Na{n, i) — a) (n, i) ■ a{n, i), Ni,{n, i) — 5^(n, i) ■ b{n, i), 

where at -a = J2l=i aia^- The left and right Casimirs are: El- El = ELi ^1^1 = ^a/2 {Na/2 + 1) , 
Eb^ • Efi ~ N\j/2 {N\j/2 +1). Therefore, the kincmatical constraints ^ in terms of the prepotentials 
mean that on every link (n, i), the number of left oscillators is equal to the number of right oscillator: 



Na{n, i) = Nbin + i,i) = N{n, i) 



(13) 



as shown in Figure ((lb)- Under SU(2) gauge transformations ^ at site n, the defining equations pT|) 
imply that the prepotentials transform as SU(2) fundamental doublets: 



aa{n, i) Aaf3{n)a0{n, i), ba{n, i) Kap{n)bp{n, i). 



(14) 



Therefore, unlike the link operators Uap{n,i) transforming at both the ends of the link by A(n) and 
A(n + i) in the prepotentials {aai'n,i),ba{n,i)) transform only at a single end by A(n). In the 
next section this simple fact will enable us to define SU(2) invariant Hilbert spaces 1-L^^'^'^\n) locally 
at every lattice site n. 

In terms of the prepotential operators we have additional U(l) gauge invariance as their defining 
equations ifTTj) are invariant under: 



ijj(n, i) A(n, i) ajj(ri, i); bl^{n + ^ ^^(n, i) bl^{n + 



(15) 
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Figure 3: The 2d prepotential SU(2) doublets a^n, i], i=:l,2,...,2d around every lattice site n shown in d 
= 2 by their Young tableau boxes □. They all transform as doublets under SU(2) gauge transformation 
at site n. 

In (fTS]). X{n,i) — expi9(n,i) is the abelian phase factor on the link {n,i). Note that the abelian 
transformations psp leave the Hilbert space of the rigid rotators, satisfying the constraints (fT^ . 
unchanged. The creation operators al^in, i) and bj^ln+i, i) create and absorb unit abelian flux at lattice 
sites n and n+i on the link (n,i) respectively. Thus the operators a]j(n, i)bj^{n + i, i), aa{n, i)bi3{n + i, i) 
create and annihilate an abelian flux line on the link (n, i) respectively. To analyze nonabelian SU(2) 
fluxes in terms of the prepotential operators, we consider the following commutation relations: 

The above relations imply that a]^{n^ i)bjj{n, i) and aa{n, i)bp[n, i) on a link (n,i), like the the link oper- 
ators Uai3{n, i), change the angular momentum j{n, i) by ±i. As an example the prepotential operators 
acting on the vacuum with j — Ocreatej = i state: {El ■ El) [aj^|0)] = | [aj,|0)] , (i?^ • Er) [b'llO)] ~ 
|[6t,|0)]. These results correspond to {El ■ El) [UcfslO)] = | [C/a/3|0)] and {Er ■ Er) [UcfslO)] = 
I [t^c(/3|0)]. More precisely, using the Wigner-Eckart theorem, the link operator can also be repre- 
sented in terms of the prepotential operators [9J: 

Uai3{n,i) = F{n,i){al^{n,i)blj{n + i,i) + aa{n,i)bp{n + i,i))F{n,i). (16) 

In (fTC|l . F{n,i) = (A^(n,i) + 1) ^ is the normalization factor where N{n,i) is defined in (fO| . Note 
that the r.h.s. of (fT6|) is U(l) invariant and also has the required SU(2) gauge transformation property 
of U{n,i) given in Acting on the Hilbert space satisfying p^ . the relation is consistent with 
^ and [C/q/3, U.ys] = [J7a/3, U}^g] — 0. The relations pTjl and ^TE\i provide complete mapping from 
the original operators satisfying ^ and the constraints (jS]) to the harmonic oscillator prepotential 
operators satisfying p2|) and the constraints ([T3|) respectively. For later convenience, we define the 
following TT operation on every link: 

al{n,i) ^ aa{n,i), bl{n,i) ^ ba{n,i) (17) 

The link operator can now be written as Uaj3{n, i) = F{n, i) ( at,(n, i)b\{n + i, z) + 7r(n, i) j F{n, i). 
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2.3 Non-abelian intertwining, abelian weaving and loop states 



The advantage of the prepotential operators is that under SU(2) gauge transformations they transform 
locally as SU(2) fundamental doublets (fH]) . Therefore, the SU(2) invariant loop Hilbert space Ti^ can 
be constructed and analyzed locally in terms of Ti^'^^^'f (n) at different lattice sites n. In the next 
section this simple fact, in turn, will enable us to solve the Mandelstam constraints exactly. To 
appreciate and elaborate on these statements further, it is convenient to relabel all the prepotentials 
and the corresponding electric fields around a lattice site n. We define^: aj,[n, i], i — 1,2, ...,2d where 
aj,[7i,i] = aj,(n,i), aj,[n, d + z] = b]^{n,i),i = l,2,..,c? as shown in Figure ([3]). Thus, instead of the 
original [E1{n, i), E'j^{n, i), Uap{n, i)] description, we now have 2d SU(2) doublet prepotential operators 
around every lattice site n as shown in Figure The SU(2) gauge transformation: 

is represented in Figure Q by a single Young tableau box on the link [n, i] with i = 1, 2, ..2d. Further 
defining J°'[n,i\ — a^[n,i]^a[n,i], we get: 

r[n,i]^ El{n,i); J°[n, d + = £;^(n, i); i^l,2,...,d. (18) 

Therefore, all possible SU(2) invariant operators at site n are given by "intertwining" (anti-symmetrizing) 
any two different prepotential SU(2) doublets: 

Lij{n) = ea(i alXn,i] a^f^[n,j] = a^[n,i]- d)[n,j], i, j = 1,2, ...,2d; i < j. (19) 

In (fT^. eafs is completely antisymmetric tensor (en = 622 = 0,£i2 = —£21 = 1) and aj, = ea^a^j- The 
d(2d — 1) intertwining operators Lij{n) in correspond to putting the Young boxes in [i] and [j] 
directions vertically to construct SU(2) singlets. In Figure (jl]), we graphically show this by joining 
these boxes with thick lines. These intertwining operators along with aT[7i,i] • a[7i,j] are the basic 
SU(2) invariant operators. Any gauge invariant operator can be analyzed locally in terms of thenj^. 
Further, the d{2d — 1) operators Lij{n) at lattice site n play the role of "SU(2) gauge invariant bricks" 
which acting on the strong coupling vacuum create the complete SU(2) gauge invariant Hilbert space 
7^S'(7(2) at n. Note that Lij{n) — —Lji{n), La — implying self intertwining is not allowed. Thus 
a basis in 7i^^(^)(n) is given by: 



|/(n) 



ha hi ■■■■ h(2d) \ 2d 



hd-2(2d~l) hd~2{2d) i)i 



h2 ha hi ■■■■ h{2d) 

\ 2d 

n (i.,(n))'-(")|0), h,{n)eZ+. (20) 

d~2{2d) J 
hd-l(2d) 

In (PO]) . Z+ denotes the set of all non-negative integers and lij{= lji,ln = 0) are d{2d— 1) SU(2) gauge 
invariant intertwining integer quantum numbers characterizing the SU(2) gauge invariant Hilbert space 
at the site n. These SU(2) singlet states can also be graphically represented by first assigning lij{n) 
Young tableau boxes individually to the links [n,i] and [n,j] and then joining (intertwining) them 



^Unless stated explicitly, the direction indices [i], [j].. within square brackets will vary from 1 to 2d. 
®See section 3 and appendix B for an explicit analysis of Tri/piaquotte- 
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(a) (b) (c) 

Figure 4: Graphical representation of SU(2) invariant intertwining illustrated for the states |Zi ), ) 
and IZ3 ). The constraint |Zi ) = |^2 ) — 1^3 ) at lattice site n represent the Mandelstam constraint in 
Figure ^ in terms of the prepotentials. The thick lines should be compared with the corresponding 
thick lines in Figure ([T]). 



together. In Figure ^ we illustrate this graphical representation for the simple case involving the 
following three vectors m d = 2: 






(21) 



The above three states |^i ) = L12L34IO), \l2 ) = L13L24IO), IZ3 ) = Li4_L23|0) are manifestly invariant 
under SU(2) gauge transformation. More expHcitly, \l{n) ) in (^0]) satisfy the SU(2) Gauss law (fTT 
constraints: 



r(n)|Z(n) ) 



2d 

E 



J''[n,i]\l{n) 



JtotailKn) 







(22) 



as [C°(n),Ly (n)] = and C°(n)|0) = 0. Thus all the states \l{n) ) in ([^ at site n are the eigenstates 
of the total angular momentum at n with eigenvalues zero. Infact, the states \l{n) ) are also eigenstates 
of individual 2d Casimirs in 



J[n, i].J[n, t]\l{n) ) = j[n, i]{j[n, z] + l)\l{n) ), z = 1, 2, 3, 4 



where, 



2d 



2j[n,i]= ^ hkin), kk{n) = lki{n) , lik{n) € Z+ 

k^i=l 



(23) 



(24) 



We note that ((M]) is both necessary and sufficient condition on j[n,i\,i — 1,2, 3,.., 2d to get SU(2) 
singlets. 

Having solved the SU(2) Gauss law, we now focus on the abelian gauge transformations (fT^ . 
As shown in Figure ([2)d), the abelian Gauss law (fT3|) states that on any link the number of a type 
oscillators at the left end is equal to the number of b type oscillators at the right end. This can 
again be easily satisfied by putting Na{n, i) — Nb(n, i) — N{n, i) abelian flux lines on every link (n,i). 
Therefore, geometrically, the SU{2) ® U{1) Gauss law demands the continuity of these U(l) flux lines 
on the entire lattice through SU(2) intertwining at every lattice site. In other words, all possible SU(2) 
invariant intertwining within Ti^^*-^^ (n) and U(l) weaving of the neighboring Ti^'^^^-' (n) is geometrically 
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equivalent to drawing all possible loops on the lattice leading to the loop Hilbert space 7i^. Given a 
configuration of closed loops on a lattice one can read off all the intertwining quantum numbers hj (n) 
at site n by simply counting the number of loop lines going from [i]*^ to [j]*'* direction. The number of 
loops passing through a link (n, i) = Na{n, i) — 2j{n, i) — Ni,{n + i, i) where 2j{n, i) is given in (|24p . 
The reverse is also true: given lij{n) V n, which are consistent with the U(l) gauge invariance, one can 
always draw corresponding closed loops. We now review the Mandelstam constraints in the basis (PU]) 
before solving them explicitly in terms of the prepotential operators. 

2.3.1 Mandelstam constraints revisited 

The basis \l{n) ) in ([20]) provides the local description of the Wilson loop basis in terms of the inter- 
twining operators and intertwining quantum numbers at lattice site n. To see this explicitly, one writes 
Wilson loop operator Tr corresponding to a loop 7 in terms of the prepotential operators using 
(fTB)) . It is clear that acting on the strong coupling vacuum they will produce states of the form |/(n) ) 
at every lattice site traversed by the loop 7. Thus, following the notation of section (2.1), the the 
Wilson loop basis can be locally characterized by Ai jn'^ — (Number of intertwining quantum numbers 
per site — Number of U(l) constraints per site) = d{2d — 1) — d = 2d[d — 1) integers per lattice site. 
This establishes that the basis \l{n) ) describing the Wilson loop basis locally in terms of intertwining 
quantum numbers is over-complete. To illustrate this further and more explicitly, we again consider 
the vectors \l\ ), |^2 ), 1^3 ) in ((2T|) and shown in Figure (jl]). Using the identity: 

(a^ • 6t)(ct . Jt) ^ (at . 5t)(;,t . Jt) _ (at . Jt)(5t . gt)^ (25) 
we find the vectors in (|2ip are linearly dependent: 

\h ) = \h ) ~ \h ) (26) 

Infact, the SU(2) identity ^ involving link operators corresponds to the identity (|25|) and the Man- 
delstam constraint ^ is the constraint (|26p written in terms of the relevant prepotential operators 
at site n. This can again be seen by writing the Mandelstam identity ^ in terms of prepotentials 
using (fTH]) or comparing Figure H^), ([TJd) and ^) with Figure (0^), (jUa) and ([It) respectively. In the 
prepotential language, the identities of the form (|25p and their various powers will make the |/ ) basis 
in (|20p linearly dependent. Infact, this is what was done in (O to get all the loop states on the two 
plaquettes related by Mandelstam constraints. Thus, at this stage, the problem of over-completeness 
of the Wilson loop basis can be analyzed locally in terms of the prepotentials. 

2.3.2 The solutions 

To solve the Mandelstam constraints we now need to focus only on a single lattice sit^ n. As discussed 
in section (2.2), the initial complete set of commuting observables at n consists of 4d angular momentum 

''Prom now onwards we will be working locally at a given lattice site. Therefore, we will ignore the site index unless 
necessary. Also, we use the notation: = X]a=l "^"I"' *]"^" * ~ 1,2..., 2d, the eigenvalues of are denoted by 
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operators: Jf , J^,i — 1, 2, .., 2d. Instead, we can also consider an equivalent basis where the following 
4d mutually commuting operators [111 I12j are diagonal: 

V 2 



csco = 



Jl-,---J2d'AJl+ 'h) -.{Jl + 'h+'h) ■■■,{Jl+ J2 + 'h J(2d-1)) , J total ■• J total (27) 



where Jl.^i - (Ji + J| + J3 + J2df and J^.^^ = (jr' + J^' + + ^2^)- The SU(2) 

Gauss law (|22|) demands J^ot^; = >^t^ota; — 0. Therefore, we drop the last two total angular momentum 
operators from the list ([?7|) . For later analysis, it is convenient to divide the remaining (4c? — 2) 
operators in CSCO (|27p into two parts: 

CSCO{I) = [Jl Jl .., J^J , CSCO{II) = [( Ji + J2)' , (Ji + J2 + J3 J(2d-i))' - J'd] (28) 

The CSCO (I) contains 2d angular momentum Casimir operators along the 2d directions and the 
CSCO(II) contains the remaining (2d — 2) Casimirs in the above chosen angular momentum addition 
scheme. On H.^^'^'^^n) the last two operators in CSCO(II) are equal because of the SU(2) Gauss law 
(|22p. We can, therefore, denote the corresponding SU(2) gauge invariant orthonormal eigenvectors by 

|jl,i2, ■•j2d; jl2, jl23, jl2..(2d-l) = J2d) = | jl , ^2 , jl2 , ^3 , il23 , , j2<i-l , jl2. . (2d-l) = j2d) ■ (29) 

The states in (|29p are characterized by maximum possible [Ad — 3) "good quantum numbers" which 
can be simultaneously measured at every lattice site. The SU(2) gauge invariant states |/ ) in (I^U]) 
are already eigenstates of of CSCO(I) with eigenvalues 2ji = Y^k=i ^ik (see (^5])). Therefore, we can 
relabel them in terms of their angular momenta: 

10-^ b'l, j2, ■■■■,j2d,jtotal = mtotal = 0) (30) 

We note that the mapping (|30p is many to one or degenerate because of the following discrete symme- 
tries of the angular momenta ji (i=l,2,...,2d) in 



In (j3ip *! 7^ 12 7^ *3 7^ H are any four different directions and r,s and t can take all possible ± integer 
values such that > 0. The advantage of the mapping ((50)) is that the states in pUj) with different I 
which are related by (|3ip are precisely the states related by the the Mandelstam constraints ([^5]) . In 
other words, the states related by Mandelstam constraints are degenerate with respect to CSCO(I). The 
reverse is also true: the degenerate states with respect to CSCO (I) are all related by the Mandelstam 
constraints (PS]) . As an example, we again consider the three loop states in (PT|) : |/i ),\l2 ) and I/3 ) in 
d=2 which are related by the Mandelstam constraints (|25p . In terms of the angular momenta: 

1^2 ) } ^ = ji2 = jis = 3i4 = 7;Jtotal = mtotal = 0). 

1^3) 

Also, given ji = ,72 = js = Ja = j, the above three possible partitions {^i}, {^2} and {Z3} given in (PT|) . 
are mutually related by (pij) . Therefore, we can lift the degeneracy and solve the Mandelstam con- 
straints by demanding that the CSCO (I) degenerate eigenbasis ([20]) to be the eigenstates of CSCO (II) 
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as well. We conclude that a complete orthonormal loop basis in d dimension is locally characterized by 
(4(i— 3) angular momentum quantum numbers and are given in (|29p . Note that the {2d— 3) eigenvalues 
of CSCO(II) are not free and have to satisfy the triangular constraints: 

\ji2..(k~i) - jk\ < ji2..k < ji2..(k-i) + jk, fc = 2, 3, (2d- 1) (32) 

along with ji2..(2d-i) ^ hd- 

2.3.3 Solving the triangular constraints 

The above characterization of the physical Hilbert space in terms of the angular momentum quantum 
numbers has been given in the context of duality transformation in lattice gauge theories [TT] leading 
to a description in terms of triangulated surfaces 12J. In this section, we further solve the triangular 
constraints, geometrically representing triangulated 2 dimensional surfaces, in terms of the intertwining 
quantum numbers {1} which geometrically represent one dimensional loops. 

In terms of the prepotentials (see Figure ([3])) the SU(2) invariant states \ji,j2, ■■■,j2d) in the map- 
ping ((5(1)) represent ^]^hk{= 2ji) Young tableau boxes on the [i]*'' hnk. Let us first consider (12) 
plane. To get the state \ji, j2i j2d', ju) from the degenerate state \ji, j2, j2d), we need to inter- 
twine (antisymmetrize) I12 boxes from 2ji boxes with I12 boxes from 2j2 boxes so that we are left with 
2ji2 boxes in the (12) plane. Therefore, 2ji2 = 2ji + 2j2 — 2li2. This process is sequential and can be 
repeated to get the eigenvalues of the CSCO(II) also in terms of the linking numbers: 

^12 = jl + j2 - jl2 
hs + ^23 = jl2 + js — jl23 
^14 + I24 + ^34 = il23 + j4 — jl234 



^l(2d) + ^2(2d) + ••■ + ^l(2rf) = jl2...(2d-l)(— j2d) + j2d - jl2...(2<i)(— 0) — 2j2<i (33) 

Given ji2 at a lattice site, the top equation fixes I12, the next line fixes ^13 + I23 in terms of J12 and 

ji23, , the last equation in ([55]) is an identity as it is already contained in ([M]) . In appendix 

A, a detailed technical calculation, involving properties of SU(2) coherent states, shows that the fi- 
nal orthonormal and manifestly SU(2) gauge invariant loop states \LS)n spanning the Hilbert space 
n^^^^\n) are: 

\LS)n = |il,i2, •■j2<i; jl2,Jl23, ■•■jl2..(2d-l) = j2d) = N{j)Y^ H 7~T (^^'^ l^*^) 1*^) 

« li. 

The prime over the summation means that, the linking numbers hj are are summed over all possible 
values which are consistent with ([M)) and ([5S)l . This summations corresponds to taking appropriate 
linear combination of the degenerate loop eigenstates of CSCO(I) related by the Mandelstam constraints 
to produce an orthonormal and complete basis. The normalization constant in (j34|) is given by: 

^(il, j2,jl2)^(jl2,j3, jl23)A^(jl23, j4, jl234) N{j2d,j2d,0) (35) 
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We emphasize that this sim- 



where iV(a,6,c) = (a+b+cliy. ' {-a + b + c)l{a - b + cy.{a + b - c)! 
pie construction (|34p of the orthonormal loop basis in terms of intertwining operators and inter- 
twining linking numbers Uj in arbitrary dimensions becomes extremely involved and complicated in 
terms of the link operators Uap and the angular momentum quantum number 

As an example, we again consider the states |/i ),\l2 ) and [l^ ) and solve the Mandelstam constraint 
(|26p. Our result immediately gives us the corresponding two independent (orthonormal) states: 

|jl = J2 = J3 = j4 = ^ , jl2 = 0) = i , 

|J1-J2=J3=J4 = ^,J12 = 1> = ^[ifs) + (36) 

It can be explicitly checked that the above two states are individually normalized and mutually or- 
thogonal. 

It is perhaps worth going back and also solve the Mandelstam constraints for the example given in 
section 2.1 with arbitrary Na and Nb in 2 dimension. In this special case: 2ji(n) = 2j2(n) = Na and 
2i3("-) = 2j4(n) = Nb- The single state shown in the top r.h.s. of (O correspond to I12 = Na and hence 
ji2 = ji + j2 — I12 = at site n. The two states in the second top line of (O represent li2{n,) — Na — 1 
and therefore correspond to ji2{'n) = 1, the three states in the next line have li2{n) = Na — 2 or 
equivalently ji2{'n) — 2 and so on so forth. The (A'min -I- 1) states in the last line of ^ correspond 
to li2{n) = if iV^ < Nb and li2{n) = Nb if Nb < Na- Therefore, these states at n correspond to 
ji2(f^) = 2ji if {Na < Nb) and ji2 = 2(ji — ^2) if (Nb < Na)- Using the relative weights in ([M]) . one 
can easily construct an orthonormal basis in 7i^^(^)(n). 

As mentioned before, the SU{2) (g) U{1) invariant orthonormal loop states on the lattice can be 
obtained by drawing all possible loops on the lattice and computing the Ad—S "good quantum numbers" 
in ((^^ and constructing ((M)) at every lattice site. All possible orthonormal loop states on lattice can 
be formally written as: 

|i'S')iatticc = I Y\. J d(l}{n,i)exp i (j}{n,i){j[n,i]~ j[n + i,d + i])\ '^\LS)n (37) 

The auxiliary angular fields 0(n, i),0 < (t>{n,i) < 2tt in ([57|l on every link implement the U(l) Gauss 
law (fT3|). Note that the U(l) Gauss law provides d relations per lattice site. Therefore, without any 
over counting, there are Af/n^ = Ad — 3 — d = 3{d — 1) "good quantum numbers" associated with the 
loop states per lattice site in the final construction ([37]). This was the desired number we started with 
right in the beginning (see equation ([T])). 



^The states \ji, j2, ji2, js, ji23, ■■■■,j2d~i, ji2..(2d-i) , jtotal,'nitotal) can be obtained by using Clebsch- 

<il2.-2d=0"ll2..2d=0 
Jl2..(2d-1) '"I2..(2£i-1) J2d"^2d 



Gordan coefficients: \ji,j2,ji2,j3,ji23,---ji2..(2d-l) = hd) = Em C^^^-^'' °'"i2..2d 



^ii2"n2^/3™3*^ii"M^/2™2 ^iii™*)- However, such approaches lead to rapid prohferation of gauge non-invariant 

Clebsch Gordan coefficients |16| forcing one to use graphical methods to avoid this problem. In contrast, the construction 
1134 I I in terms of gauge invariant intertwining numbers (not angular momentum) is simple and bypasses this problem. 
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2.4 Matter, loops and strings 



The inclusion of fundamental matter fields is also natural in the prepotential formulation. It simply 
increases the number of intertwining operators at every lattice site. For simplicity, we introduce 
scalar matter field operators: {(j)a{n),(j)*^{n)) and their conjugate momenta (71^(71), 7r*(n)) respectively. 
They are neutral under U(l) gauge transformations and transform as doublets under SU(2) gauge 
transformations: 0Q(n) AQ,^(n)0;3(n), tt* (n) ^ AQ,^(n)7r^(n). The canonical commutation relations 
are: 



['^a{n),(j)l3{m)] = [7r*(n),0* (n)] = - i SapSnm: 

[■Ka{n),TTi3{m)] = [(j)a{n),cj)i3{m)] = 0, 
[7ra(n),7r^(m)] = 0^(m)] = 0. 

The matter creation and annihilation operators are defined as: 



Like the prepotential operators they satisfy: 

aa(n),at(TO) = 1, 
Za{n),bp{m) = 0, 



1, 
0. 



Under SU(2): 



(38) 



(39) 



(40) 



(41) 



they transform exactly like prepotentials p4p thus putting matter and gauge sector on the same footing 
under non-abelian gauge transformations. Therefore, we now have to construct SU(2) singlets out of 
{2d + 2) types of SU(2) doublets of creation operators per lattice site instead of 2d types as in the case 
of pure gauge theory. The orthonormal SU(2) gauge invariant states are now characterized by: 



|jl,i2,il2, ja, jl23, il2...(2rf-l)i j2d, jl2...2(i, j(2d+l)iil2...(2rf+l) — j(2rf+2))- 



(42) 



Further, the iterative method in appendix A again goes through and the states ([42|) can be easily 
constructed in terms of intertwining operators which will now involve matter creation operators (|39p 
also. 

As the matter fields transform like SU(2) doublets and are neutral under U(l) gauge transforma- 
tions, they provide end points for the abelian fiux lines leading to a gauge invariant description of 
lattice gauge theories in terms of loops and strings. Note that including scalar matter in the funda- 
mental representation amounts to d ^ d + 1. The extra quantum numbers required per lattice site 
due to the inclusion of matter is [4 (d + 1) — 3] — [4(i — 3] =4, as expected. 
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3 The dynamical issues 



To discuss dynamics of loops, we consider pure SU(2) lattice gauge theory Hamiltonian [5]: 

i/ = ^^i?2(n,j) + l^(2-TrC/n) (43) 

where, 

Ua - C/piaquctto = U(n,i)U{n + iJ)U''{n + j,i)U''{n,j), (44) 

g is the coupling constant and E^{n, i) = E\{n, i) = Eff {n + i, i) (see (|S])). The loop states discussed in 
section (2) trivially diagonalize the electric field term with eigenvalues "^i^^pg j{n, i) {j{n, i) + 1) where 
Sloops denotes summation over all the links on the loops. The electric field term in this loop basis 
is like potential energy term and counts the number of abelian flux lines on that link. The plaquette 
term acts like kinetic energy term: it makes the loops fluctuate over the corresponding plaquette by 
creating and destroying the abelian flux lines (fTH)) . 



3.1 The loop dynamics 

As mentioned in the introduction, we restrict ourselves to d=2 and generalize the results to arbitrary 
d dimension at the end. We consider a plaquette abed as shown in the Figure ^ with the four edges 
(ab), (be), (cd), (da) denoted by /i,Z2,^3,^4 respectively. Using and P?)) . we write the gauge 
invariant plaquette operator over abed in terms of the prepotentials: 



T^rUabcd = Fa, 



bed 



(at[l] . at[2])^ (at[2] • 5t[3]), (at[3] • ~a^[A])^ (at[4] • ~a^[l])^ + ^ nik) 



ij>i—l ij>i,k>j — l a/37(5— it 



where Fated = F{li)F{l2)F{l3)F{l4). In there are sixteen SU{2)($U{1) gauge invariant terms 
which are produced by substituting (|16p in (|44p . The first plaquette operator on the r.h.s. of the 
top equation in (j45p is written explicitly in terms of the prepotential intertwining operators at lattice 



a^[2] 



a + 

a ni 



a"^[3] 



aVi 



a^[2] 



+ b 

a [31 



Figure 5: The plaquette abed with the corresponding prepotential operators from Figure ([3]). 
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Figure 6: From loop kinematics to loop dynamics, (a) The angular momenta satisfying SU(2) Gauss 
law on the dual lattice, (b) The matrix elements {jabcdlTrUabcdljabcd) with X = ^ (also see [12]). The 
six tetrahedron are the six 6j symbols in (|48p . The unchanged j lines represent the delta function and 
a • over a j line represents the factor (— l)-'n(j) leading to Mabcd in (05]). 



sites a, b, c and d in Figure (|S]). It increases angular momenta by ^ on the four links of abed and 
therefore we represent it by The rest of the 15 plaquette operators in (|i5]l are generated 

by the action of tt on it. The single tt operation at any of the 4 links (li, I2, h, li) acting on 

produces 4 terms: H ^^^,H^ ^^,H^^ ^,H^^^ . Similarly, the double tt operation produces 

6 operators: H ^+;H ^ ^_,H h'-^H i ' -^+H • There are four 3 tt terms in (H5)) : 

,H I ,H ^ ,H ^. Finally, the 4 tt operation on all the four links produces a single 

term: H . Note that (^H^^^^^ = H_^_p The advantage of the form (j45]) is that now this 

magnetic field term can be analyzed locally at a, b, c and d in the manifestly SU(2) gauge invariant 
way. To compute the action of Uabcd, it is convenient to label the loop state (j34p at abed as \ jabcd)- 

^ IfJ^j^fM ® \jlj',jWl2) ® |j?J2J3J4Jr2) ® IJI J2 J3i4 ■ (46) 

The U(l) Gauss law demands: 

jI ^ Ji , - Jl EE n , J3 = = J3 , 3i = J2 ^ J4 . (47) 

The matrix elements of TrUabcd are computed directly in the loop basis (I46p in appendix B using 
generalized Wigner-Eckart theorem and Biedenharn-EUiot identities. The final result is: 



{jabcd\TrUabcd\jabcd) — 'Mabcd 




b :;b 1 

2 



Jl2 Jl2 



.n j1 



■d -d 1 
J12 J12 2 



(48) 
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In (gHl), Mabcd = DabcdNabcdPabcd where: 

Nabcd = n (ji , ji , j2 , J2 , ja , js , j4 , ji , jl2 , jl2 i Jl2 : Jl2 ) (49) 

In ()49p . Dabcd describes the trivial (5 functions over the angular momenta which do not change under the 
action of the plaquette operator (|45p . Nabcd and Pabcd give the corresponding numerical and the phase 
factors respectively. The multiplicity factors are: Il{x,y,...) = y/ (2x + l){2y + 1)... and A{x,y,z) 
represent the phase factors associated with a triangle with sides x, y, z: A{x,y,z) = (— =^ 
A{x,y,z) = ±1. The matrix elements describe the dynamics in the loop basis and can be 
geometrically represented by the Figure ^jp). This dynamics contains three physical discrete angular 
momentum loop co-ordinates numbers per lattice site. The matrix elements (j48p have been obtainecl^ in 
the context of dual description [T^El] of (2 + 1) dimension lattice gauge theory in terms of triangulated 
surfaces. It is easy to see that the matrix elements in (1^5)) are symmetric: the 6j symbols satisfy 















X u 1 


\-\ 




X U ] 




y V J 




I y 


y V j 



and the factors Dabcd, Nabcd and Pabcd are individually symmetric 

under jabcd jabcd- The matrix elements in (j48p are also real as TrC/piaqucttc is a Hermitian operator. 
This reality can again be easily seen as the 6j symbols, Dabcd, Nabcd, A{abc) are themselves real. The 
remaining two phase factors (— l)Ji+J2+ii+i4 g^^^j 1^3+^4-1-^2+^3 in Pabcd are real because {ji,j2,ji, jl) 
and {ia, ji, j^, jf) are the coordinates of the loop states at b and d respectively. Therefore, implies: 
ji + 32+ ii+ ii = Integer, jz+ji+ji + ii = Integer =^ Pabcd = ±1- At this stage, before generalizing 
the loop dynamics to arbitrary d dimension, we cross check the d = 2 result . As the six 6j symbols 
and the 5 functions in Dabcd are geometrical in origin, we only need to check the numerical and the 
phase factors Nabcd, Pabcd respectively. For this purpose, we replace TrUabcd by the identity operator I. 
The computations in appendix B imply that now we only have to replace ^ in each of the six 6j symbols 

in gHI) and in Pabcd in (011) by 0. Using the value | " ^ | {-^T^^'^'^Sa-aK.b (n(a,6))"\ we 
get: 






Jl Jl J4 J I h k jl J _ ^jlJl^j2j2^j3j3^34jJjl2rjlJjt2Jt2 



J12 J12 M 1 Nabcd Pabcd 

h js ji jl ii ii ji J 



(50) 



Geometrically, the equation (|5(I)l corresponds to putting X = in the Figure ([Sb). It implies 
{jabcd\^\ jabcd) = ^JabcdJabcd Confirming the numerical and the phase factors in (j48p . We now write 
(|48[) in a more compact form which can be directly generalized to higher dimension. Henceforth, we 



^Our phase factors in I l49t are different resulting in real and symmetric matrix (jabcdl'^^Uabcdl jabcd) in II48I I. 
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Figure 7: The 18j ribbon diagram representing exact SU(2) loop dynamics without any spurious 
gauge or loop degrees of freedom in d=2. The interior (exterior) edge carries the initial (final) angular 
momenta and the six bridges carry the angular momenta which are invariant under the action of Trlla ■ 
The six bridges along with the respective four angular momenta attached represent the six 6j symbols 
appearing in (|51[) with ji = jizL ^. 

ignore Dabcd representing trivial 5 functions in (j48p . We write: 

X y J4 J4 ji2 J y 3l2 jl2 h J 



Jl2 Jl2 ^ 

h h ii 



J3 J3 X 
h 32 ju 



J2 32 X 

'lb -b -b 
3l2 3l2 Jl 




bed 



31 34: 

A a 
3l2 

31 k 



Ad 

33 



Ad 

Jl2 



-Ad 
3l2 



Ad 
32 



33 32 
3l2 

33 32 



3i 



3l2 



3l2 



18j coefficient of the second kind 



1^ 
J4 



i=l 



(51) 



The 18j symbols in ([STjl arc shown in ([7]). Note that Pahcd (= (^1)*^^^) in ((49|) is precisely the phase 
factor needed to define ISj symbol [26 in ([5T|l . Further, the 12 triangular constraints in ([5T|) at the 
12 vertices of the ribbon diagram in Figure ([7]) are already solved in terms of the linking numbers. 
Therefore it is only the value of the 'inj = 18j (n = 6) symbol which is important. The form ([5T|) also 
makes reality and symmetry of (jabcd |TrUabcd|ja6cd) manifest as 3nj symbols of second kind are real 
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Figure 8: (a) The plaquette abed in the (/, K) plane in d dimension. We choose I < K, 1 < I < d 
and 1 < K < d, (h) The angular momentum addition scheme at site x (=a,b,c,d). Note that jf and 
appear twice in the scheme. The S functions are subtracted in (|53p to avoid this double counting. 
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Before going to arbitrary dimension, we make the following simple observation. Let AiV^;, x — 
a, b, c, d denote the number of angular momenta appearing in the loop states \jabcd) m (j46|) which 
change under the action of the plaquette operator TrUabcd at lattice site x. In the present, d = 2, case: 



AN,, 



2, UIj^); an, ^3, (j'^jIjI^); AN, 



3; {3 1^ it 13 12) ■ 



2, [jIjD- AN, 

The U{1) identification (|T7|) implies double counting on each of the 4 links of the plaquette abed. 
Therefore, the number of angular momenta which change under the action of the plaquette in the (12) 
plane: AiV(12) = ANa + AN, + AN^ + AN, - 4= 10 - 4 = 6 = n. This analysis will be useful to 
generalize the loop dynamics to arbitrary dimensions below. 



3.2 d dimension 

It is clear from the previous section that the loop dynamics in d dimension is also given in terms of 3nj 
symbols. However, in arbitrary d dimension, n will depend on the orientation of the plaquette. We 
will now compute n. We consider the plaquette abed in the {I,K),I < K plane as shown in Figure 
Like in d = 2, we consider the loop states over the plaquette abed: 

\jabcd) = \LS)a ® \LS)b ® \LS)c ® \LS)d (52) 

where \LS)x=a,b,c,d = lif jjf 1^127^3 'il23' ■■■■^ 3l : 3l2...I^ 3kt 3l2..I..K^ ■••■'.?2d-l'.?'l2..(2d-l)(~ 32d)^32d^ ^) 

We now have to count the the number of angular momenta in (I52p which change under the action of 
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Figure 9: The 30j ribbon diagrams representing exact SU(2) loop dynamics without any spurious 
gauge or loop degrees of freedom in d=3: a) (12) plane, b) (13) plane, c) (23) plane. The angular 
momenta Ji,j2, js and j4 are as shown in Figure ([5^). The initial (inner) and final (outer) angular 
momenta differ by ^. 

the plaquette operator Trllabcd in the (IK) plane. With the choice I < K {1 < I < d,l < K < d) , we 
have: 

ANa^2+(K - I)-di,i, ANb^2+{d + I)-K, 

ANc^2 + {K-I)- SK,d: ANd = 2+{d + K)-I - 5is - Ska (53) 

This implies: 

AN{IK) = ANa + ANb + AN^ + AN^ - A = 2[2 + d + (K - I) - (5/,i - dK,d] = n{IK). (54) 

Like in d =2 case, we have subtracted 4 in fM)) because of U(l) gauge invariance. Note that for d=2, 
AiV(12) = 6 and for d=3, AiV(12) = AN{13) = AiV(23) = 10. The d = 3 loop dynamics is explicitly 
shown in Figure ^ where we have used the notations from Figure ([Hi), i.e.: 

{I^1,K^2) = (12) plane : = jI = n, - j| = j2, 3% = Ji = .?3, = 3l = J4, 
(/ = 1, if = 3) = (13) plane : = jl = ji, j| = jl = 32. 3% = Ji = 3z.jt = H = .H, 
(/ = 2,K = 3) = (23) plane : = jI = ji, jI = = j2, j| = = js,.!^ = it = Ja- 

It is clear from (|54|) that in higher {d > 3) dimension AN{IK) depends on the orientation of the 
plaquette. The corresponding 3n(/, K)i symbol describing the loop dynamics in the above angular 
momentum addition scheme can be easily written down. 

4 SU(N) prepotentials 

We now briefly discuss the extension of the ideas in this paper to the SU(N) gauge group. The SU(N) 
generalization of the SU(2) Jordan-Schwinger mapping (|lip has been done in 18J in the context of 
SU(N) coherent states. We define the left and right SU(N) electric fields through the SU(N) prepo- 
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tentials: 

N-l N-1 .(jr , N-1 N-1 , 

El^J2 - E ^Hr]^a[r]; E% ^ E^r] - ^ bHr]^b[r]. (55) 

r— 1 r — 1 r—1 r—1 

In ([55|) . the index r(= 1, 2, (A^ — 1)) covers all the {N — 1) fundamental representations of SU(N) 
group and all the link indices (n,i) are suppressed. Again, the {N — 1) SU(N) Casimirs on the left and 
the right ends of any link (n,i) are the prepotential number operators: 

Na[r]=a^[r]-a[r], Nb[r] ^ [r] ■ b[r]:, r = 1,2, .., {N - 1). (56) 

The SU(N) kinematical constraints involving SU(N) Casimirs ^\ are now in terms of the prepotential 
number operators: 

Na[r]in,i)^Nb[r]{n + i,i), r ^ 1,2, ..{N - I). (57) 

The (A^— 1) constraints (|57|) are analogous to the single constraint (fT3|) for SU(2). The defining relations 
(|55p for the prepotentials and the constraints (j57p imply the following abclian gauge invariance: 

a^[r]{n,i) ^ expi9[r]{n,i) a^[r]{n,i); b'^[r]{n + ^ exp — i9[r]{n,i) b^[r]{n + (58) 

Thus the SU(N) prepotential formulation will have SU{N) (g) (C/(l))^ ^ gauge invariance leading to 
(A^ — 1) varieties of loops. Therefore, using the ideas of SU(N) Schwinger bosons and the corresponding 
coherent states [H], like in appendix A for SU(2), it should be possible to find an orthonormal loop 
basis and corresponding loop dynamics for SU(N) lattice gauge theory as well. The work in this 
direction is in progress. 

5 Summary and discussion 

In this work we developed ideas and techniques to formulate Hamiltonian lattice gauge theories exactly 
and most economically in terms of loop and string degrees of freedom. This required systematically 
solving Gauss law constraints, Mandelstam constraints and finally triangular constraints. The ap- 
parently highly non-local and formidable Mandelstam constraints in terms of the link operators were 
cast and then solved locally in terms of the prepotential intertwining operators. Infact, one of the 
motivations for this work was to develop manifestly SU(2) gauge invariant techniques involving gauge 
invariant local intertwining prepotential operators and intertwining/linking quantum numbers having 
direct geometrical interpretation in terms of loops. The resulting simplifications have been emphasized 
in the text. Note that the loop states as well as loop dynamics were constructed directly in terms of 
the above, without using Clebsch-Gordan coefficients which are not gauge invariant and also do not 
have direct interpretation in terms of loops. The final loop dynamics, i.e., the matrix elements of 
the magnetic field terms in between the two loop basis states, are found to be real and symmetric 
and are given by 3nj symbol in arbitrary dimension. Therefore, this loop space description of gauge 
theories is also a non-abelian dual description [l^] where the effect of compactness of the gauge group 
is contained in the integer intertwining or half-integer angular momentum quantum numbers labeling 
the loop states. In the simpler context of compact (2+1) and (3+1) U(l) gauge theories such duality 
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transformations are known to isolate the topofogical magnetic monopole degrees of freedom leading to 
confinement |27| . It will also be interesting to develop a systematic weak coupling {g — s- 0) loop per- 
turbation theory near the continuum. This perturbation theory should encapsulate the global gauge 
invariant loop fluctuations as opposed to the local gauge field fluctuations which is the case with the 
standard perturbation theory. The issue of color confinement and vacuum structure will be of special 
interest. The work in this direction is in progress and will be reported elsewhere. The present prepo- 
tential approach has also been found useful to study spin networks and analyze the spectrum of the 
volume operators in lattice loop quantum gravity [55]. Finally, the maximally reduced loop basis and 
the corresponding matrix elements should also be useful for numerical diagonalization. 

This work is dedicated to the memory of late Prof. C. K. Majumdar. 



Appendix A 

In this section we explicitly construct all possible orthonormal loop states in terms of prepotentials 
intertwining operators. We appropriately interpret, modify and generalize the techniques developed 
in [17| to add angular momenta in terms of Schwinger bosons or equivalently prepotentials in our 
formulation. The basic idea is that angular momenta can be combined directly in terms of prepotentials 
and the d{2d~ 1) intertwining operators (|19p by taking certain direct products of SU(2) coherent states. 
We explain this idea in d = 2. It's generalization to arbitrary dimension is then obvious and done 
next. The SU(2) group manifold is characterized by a doublet of complex numbers (2:1, Z2) with the 
constraint: |zip + |z2p = 1. The SU(2) coherent state in the spin j representation are given by [TH]: 

( \J ~j~ ^Tl / \j — /?l 

\z), = n ^ T 1^'"^) ^ E <^^»(^) IJ''^) (59) 

— ^ . V (7 + to)!(? — ml! 

m——j V v-' ' / v-' / m——j 

In terms of Schwinger bosons, 

\j,m)^ )y y ^, 10,0). (60) 
The generating function of SU(2) coherent state is: 



E 

3 



<i>j{S)\z)j = exp {S z ■ a^) |0) = exp (s {zia\ + Z24)) |0) (61) 



where ^j{5) — . The states |j, m) can be extracted by comparing the terms with coefhcients 5"^^ 
on both sides of (gl]). To add the 2 angular momenta, J'^[l] = at[l]^a[l] and J°[2] = at[2]^a[2], 
we consider direct product of the generating functions of two SU(2) coherent states defined over the 
complex planes {xi,X2) and (yi,?/2) respectively. 

k> ® \y) = E ® \y)32 = exp {x -a^lM+y a^2]) |0) (62) 

We apply the differential operator involving a triplet of complex parameters {Si, (^2, (S3) and a complex 
doublet z{= (zi,Z2)): 

exp (^63 {dx ■ By) + 5i (z • dx) + 62 {z ■ dy)j (63) 
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on the both sides of (|62p and put x = y = to get |17j : 

E Wk>^{f = exp(<53 a^[l] ■ a^[2] + z ■ at[12]) = expf^S^ + z ■ at[12]) |0) (64) 

jlj2jl2 

where L12 is the intertwining operator in the (12) plane and 

ai[l2]^S^al[l]+S2al[2] (65) 



(2jl2 + 1) 



(^^-|^)il-i2+J12 ^J2)~-?'l+-'2+j'l2 (^§^yi+h-jl2 



[{ji ~ h + ji2)!(-ii + .72 + ji2)!(ji + h - ii2)!] 



(66) 



In Iz)^!^^ is the coherent state in the ji2 representation of the combined angular momentum 



U1J2 
I hi 

(/[!] + J[2]),i.e: 

+jl2 



\z 



^3\h — 
ii2 ~ 

'ni2 = -jl2 



'/'il2mi2(^)|jlj2jl2mi2) (67) 



Note that (j\ — 22 + ^12), (— ji + J2 + J12) and (j\ + ji — J12) are all non-negative integers due to 
angular momentum addition rules. The coherent state generating function (|64p in the representation 
ii2 (|ji — 32\ < ji2 < ji + ji) is the generalization of the generating function (j6ip in the case of single 
angular momentum. To illustrate this, we put 82 = St, = 0{=^ j2 = = ji2) in (|64p to recover (|6ip . 
Now to project out Jtotai — J[^] + J[2] ~ ("gauge invariant state in d = 1") state from the left hand 
side of ([64]) . we put: 



61^0, (52 = ^ ji = j2 = j, jtotai = jl2 = 0. 

Now takes the form: 



|jl = J, J2 = J, Jtotai = J 12 = 0) = Af{j) |0) (68) 



mj) = - 32 = J, J12 = 0) = -=== (69) 

V (2j + 1) 



In dM]) ?i2 = 2j and 



In ^(J1,J2,J12) = V (TTtItStTT [(-^'i +J2 +Ji2)!(ji - J2 + ji2)!(ji +J2 - Ji2)!]^ are the nor- 
malization constant so that (ji = j, j2 = j, jtotai = ji2 = 0\ji = j,j2 = j Jtotai = ji2 = 0) = 1. The 
U(l) Gauss law makes j site independent. These states are the trivial examples of the "loop states" 
in d=l. We now generalize the above simple construction in d=l to d=2. To add the third angular 
momentum corresponding to [n, 3] we rewrite (|64p with z ^ x and take the direct product with 

\y) = X] '?^J3m3(y)lj3,m3) = exp{y.a'^[3])\0) 

We now apply the operator (|63p with 5 ^ (t: 

exp (^as {dx ■ By) + ai {z ■ d^) + (T2 {z ■ dy)^ (70) 
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to get: 



where, 



E 'l>nnnAS)^n.nn.M\^)]lT''' ^ |0) (71) 



a''[123]=aia''[12]+a2a''[3] 
L123 = S3 a^[l].at[2] + CT3 at[12].at[3] = + ^3 {SiLrs + &2L23) ■ (72) 

and 

+il2 3 

I^)ii2il23 = E 0ji23mi23(^)|jli2jl2j3jl23'™123) 

^123— — J123 

is the SU(2) coherent states in the combined angular momentum J[123] — ^^J'[l] + J[2]^ + J'[3]^ 
representation with |ji2 — jsl < J123 < ji2 +^3- The operator L123 contains the intertwining operators 
in the (12), (13) and (23) planes. Note that, like in the previous case, putting cr2 — (^3 — 0(=> — 
0, ji23 = ji2) in ((7^ we recover ([M)) . Therefore, it is a sequential process. Repeating the same steps 
once again in the 4*^ direction with the prepotential operators a^[4], we finally get: 

E "^nn,.. mn.nn.. {^)'^3...,.n... {p)Vy^^:''''''''' = exp (L1234 + z.a^ [1234]) |0) (73) 



J1J2J12J3J123J1234 



where, at[1234] = piat[123] +p2a'^[4] and L1234 contains intertwining in (12), (13), (14), (23), (24) and 
(34) planes: 

ii234 = <53 a^[l]-a^[2] + fT3 at[12].at[3]+p3a^[123]-at[4] 

= (53 (-^12) + 0-3 ((5iLi3 + (52^23) + P3 (0-l<5lil4 + 0-1(52^24 + 0-2X34) (74) 

The gauge invariant states can now be projected (like in d=l case) by choosing: 

P=(0, 0,P3) => jl23=i4, jl234 = 0. 

We now compare the coefficients of ((^i, (^2, <53), (ci, (72, 173) and p3 to get the all possible manifestly 
gauge invariant orthonormal states at site n. After some algebra we geio: 

|jl, j2, jl2) = j2,jl2, (i3), jl23 = Ui)diJtotal = il234 = O^rUtotal = "li234 = 0) 

- AA(,) E'R^TTT-^IO) (75) 

'13>'14 i,i , ^ 
^23''24 

In ((751) . -^(i) is constant given in terms of N{ji, j2, ju) in ((M)) : 

A^(j) = A^(J1,J2,J12) A^(jl2,j3,j4) ^(J4,J4,0) (76) 

so that the states in ((75)1 are normalized to unity. The d=2 results ([75)1 and (|75)) are analogous to 
d=l results given in (|68p and respectively. The difference is the summation over the linking 



^Note that we have solved the SU(2) Gauss law at a site n and have ignored the site index from j" and P 
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numbers in (j75p which was missing in (|68p . This is because now there are many possibihties of hnking 
or contracting the 4 types of Young boxes representing [1] , [2] , [3] and [4] mutuaUy to produce 
jtotai = states. These contractions must be such that the states ([75]) are the eigenstates of the 
CSCO(I) as wch as CSCO(II) hstcd in In the construction ([75]) the hnking numbers have to 

satisfy: 

2d 

^Z„=2j„ i = l,2,..,2d. (77) 

and 

h2 = ji + h - ji2 
hs + hz = ji2 + is - ji23 = ji2 + is - ji (78) 

hi + hi + 134 = jl23 + j4 — jl234 = 2j4. 

ensuring that the states are common eigenstates of CSCO(I) and CSCO(II) respectively. The con- 
straints ([77)1 and ([75)1 on the hnking numbers are denoted by the / over the summation sign in ([75)1 . 
The last constraint in ([75)) represents the partition of ^4 in terms of the linking numbers and is already 
contained in ([77)) . Note that we have combined angular momenta directly in terms of kj without using 
any Clebsch-Gordon coefficients containing magnetic quantum numbers. The generalization of these 
techniques to arbitrary dimension is a sequential process and the results are given in ([M)) . 



Appendix B 

In this appendix, we compute the SU(2) loop dynamics. As mentioned before the electric field term in 
(|43p simply counts the abelian fiux lines without changing the loop states. The action of the magnetic 
field term in ([43| on the loop states is non-trivial. The calculations of the matrix elements of TiUa in 
d=2 are presented below. It is convenient to define the square root of the multiplicity factors: 



which will occur repeatedly below. 

The matrix elements of an intertwining operators L12 — (a^[l] • a^[2]) in the corresponding angular 
momentum basis j2,ii2 7 "^12) are given by the generalized Wigner Eckart theorem |251 126] : 

{ji,j2,ji2,mi2\ (at[l] ■a'^[2]) | ji, j2, ji2, m-12) = V2 {ji, j2, ji2,mi2\ (a'l'[l] (g) a^p])" |ji, j2, ji2, mi2) 
= V2(-l)&-'^-)n(ji2,Ji2,0) f I ]\ j2 32 \ ^(ji||ani]lb-i)(J2||at[2]||j2) 

\ -77112 mi2 / - 

[ J12 J12 

In ([7^ . {a^\\\ ® a^[2])g = m=±\ m°^rn}^m ^ith i+i = ^i, a_i = 02. The reduced matrix 

elements of the prepotential operators are given by: 

(j||a||j) = n(j, j) , Q\\a^\\j) - n(j, J) <5.^^.+ i. (79) 
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The coefficients ^ 
Using the values: 



31 J2 J3 

mi m2 fTi3 



jl2 ji2 
-77112 TOi2 




represent 3j and 9j symbols [25] respectively. 



(-1)- 



-Ji2+mi2 



Jl Jl 
32 32 
312 312 



n(ji2) 




we get: 



(jl , j2, jl2 , "I12 |Li2 I jl , j2, jl2 , mi2) = 



^Jl2jl2'^"n2,mi2(-lF''?7(jl, j2) 

(ji||ani]l|ji>(J2||at[2]||j2). 




(80) 



In ([SO]) , the phase factor ri{a,b) = (— 1)''+^+2. Note that one can also apply the above intertwining 



operator directly on the loop basis ([75]) to get it's matrix elements (|80l) algebraically. Infact, this has 
provided an independent check on the "master formula" (|80p which is crucial for the computations 
below. It simply states that the intertwining operator L12 — a^[l] • 5,^^[2] increases ji and j2 by i. 
Further, as L12 commutes with the ( J[l] + J[2])^, the matrix elements are diagonal in ji2 and r7ii2. 
The matrix elements of the intertwining operators in the geometrical form ((80|l also tell us that all the 
16 terms in J2ai3-yS=± ^a/SjS 651) differ only in their reduced matrix element structures. Therefore, 
we need to compute the matrix elements of only a single term in (j45p , providing enormous simplification 
at the algebraic level. In the following calculations, we choose this term to be the first term in 
associated with the plaquette abed in Figure ([5|): 



+++ 



Fa,.a (ani] • «n2]), («n2] • «n3]), (^nsi • ~a^[^), (an4] • aniD.^^a&c 



(81) 



As mentioned in the beginning of this section, we will work in d = 2. The results are then easily 
generalized to arbitrary dimensions d and are given in the equation (|5ip . In computing the loop 
dynamics below, it will be often convenient to change the angular momentum addition schemc0: 



|jl, j2, js, j4, jl2) = |jl, j2, jl2, js, jl23(= j4), j4, jtota; = j(123)(4) = 0) 
|(jl, j2)jl2, (ja, j4)j34(= jl2)Jtotal = j(12)(34) = 0) = |(jl, j2)jl2, (js, j4)jl2) 



(82) 



^^The relation I I82I I is proved by writing: 

Iil,j2, j3,i4,J12> = ^ C^^ 



0,0 

23™123>i4"l4 



/—I * J121112J 
all m js4 



, \jl , j2 , J12 , mi2 ) lis , j4 , J34 , "134 > = ] (jl , j2 )jl2, (js , i4)jl2 > • 



We have used: 
0,0 



C,, 



(-1):'4- 



■^123"^123 ■■?4 "^4 



'^■?123'-^4 '^"^123' 



(-1)^4 



n(jj) .'l2'"l2J'3™l 



qO,0 (-rJl2-'"l2 
Jl2'"l2 Jl2-'"l2 iS"^t,jS"^i' 
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The equivalent scheme on the right of ([82)1 simpUfies the algebra. We also note that the normalization 
operator Fabcd in has simple action on the loop states \jabcd) defined in ([i^ and (|T7| : 

Fabcd\jabcd) = TTTT : : -^\jabcd) (83) 

nUl,j2,j3,j4) 

Therefore, we only need to compute the matrix elements of the intertwining operators in ()8ip in the 
orthonormal loop basis given in (|75p . 

Loop dynamics at a: 

In H^^^^ above, the intertwining operator at a is {a^[l] ■ a^[2fj^. Using ([50]) . one directly gets: 

I -! 1 ! |(Jrikni]ibT>(j2ii«n2]iij2> (84) 

[ 32 J2 Jl2 J 

The intertwining operator {cl^[1\ ■ at [2])^ increases the SU(2) flux on the links li as well as Z4 of Figure 
([Kjl. Note that this information is contained only in the last two reduced matrix element terms in 



Loop dynamics at b: 

The intertwining operator at b in ([81]) is (at [2] • at[3])j^. To compute it's action at b, we write the 
loop states ([M|) in terms of the basis states which diagonalize (J[2] + J[3]) 



2. 



nb 



\Ajiiji3\^ = (-i)^^^+^^+^^+^-^En(j?2,j^3) 'I 



^23 



2z J123 ^23 



- (-i)'^^+^^+^^+^^^En(^'^23)j 'HlJ?^^j^^» (85) 

jig I J23 J 

In ((SSD, |jtj2j1iti23) = Iji'-(j2i!)'i23>im'i4':im4 = ™?234 = 0) and we have used ^123 = H- Note 
that the phase factor (—\^('h+3i+h+'h) in ([55[) is real because of the triangular constraints on the 
angular momenta or equivalently (|24p . Now we use: 



|J?J2J3J4J23) - (-1)'^'^ b2J3i4i?i23 = jtl) (86) 

and ([HOI) to get: 

(jj J2 J3 J4, J12I (an2] • «n3]), b^ j^, jt^> = (-i)^^+^'"^^"^^^(j2 Ja)'^,? j?^,tjin(jj2, J12) 

E(-)^^^^n(,t3)| t 11 'I '\ i 11 'I f |(jlll«n2]|b2^)(jll|an3]|b3^) 

t Ji Ji J12 ) y H J2 J12 J t J2 J3 2 J 

The summation over J23 in the last line above can be performed using Biedcnharn-EUiot identity [25j : 



E(-irtf(^)( 



a b X ]\ c d X ]\ e f X \^^_^yrl PQsl^lpqs] 
cdp\\e fq\]bas\ \ead\]fbc\ 

r = {a + b + c + d + e + f+ p + q + s). 
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Finally, the loop dynamics at lattice site b is given by: 



Gin Jill Jl2\ («n2] • 5t[3]), bt, Jt Jt2> - i-^y^'''^V{jlj%{j'l2Jl2)S,^J^S^iJi 



n(jt2,J?2)< 

Loop dynamics at c: 



Ab -Ab 1 

J12 Jl2 2 

Js is il 



Ab -Ab 1 
J12 Jl2 2 

^2 J2 Jl 



{j2\\aHmj2){jl\\aHmjl) (87) 



At c, we use (f82|) to write 

lili ^2; JS' ^4; J12) = I0ijJ2)Ji2; UsJt) j Ju) = l)^"*^^ I (js ' ^1)^12 ; (jli2)jl2) 

to get: 

(j?, j1, j1, J4, Jr2i (0^3] • 0^4]), bT, jI, J4^ jr2> = (-lF-77(J3^Jl)%J^ 

•(J3iian3]iij3)(j4ii«n4]ibi) 



^3 is 5 
if J4 J 12 



Loop dynamics at d: 

To compute the loop dynamics at d, we write: 

I -d -d -d -d -d \ I -d -d -d -d -d -d / -d \\ 

\JiiJ2jJ3tJ4tJ12/ — \J1tJ2tJi2tJ3tJ4U34[—J12)/ 



Jl 



= E(-l)''^-*'''"^'^n0f„J3^4,J-1^4,J2^3)<| Jt 

= (-i)^^"En0f2,jf4)(-i) 



Ad -d _ -d 
Jl4 J23 — Jl4 



Jt2 

Ad _ -d 
J34 — J12 





(88) 



I -d -d -d -d -d -d \ 
( \J4tJ1 1 Jl4: J2 1 J3 ' J23/ 



J14 



jf+ii+ii+ii J -^1 -^'2 -^12 1 I -d -d -d -d d d \ 
' ' ' M .rf Mj4>Jl:Jl4>J2:J3:J23) 
I Ji J 3 J 14 I 



(89) 



In 



Jt 


ji 


jl2 












jl2 


A4 


rL 






have used <( jf, ^ (-1)'"^'"^'"^'" (n(ji2,jY4)) 



-1 1 ji jt2 

1 Ad -d -d 
I J3 J3 Jl4 



Finally, using (|89p and the Biedenharn-Elliot identity, the dynamics at d is given by: 

GtJijiJiM («n4] • ~aHl]), , J2^ J3, J12) = -h^Y'^^'hijfJtMjt2Jt2)S,^J^S^iJi 



n(jt2,jT2 



■d -d 1 

J12 J12 2 

'Ad -d -d 

Jl Jl J2 



^-l ^'l i |(illl«n4]|b4)(jill«ni]lbi) (90) 
J4 J4 J3 J 



Loop dynamics at abed: 

We now weave or glue the dynamics at a,b,c,d with the help of U(l) Gauss law ()47p : — J3 = ji, j\ — 

jl = 32, Jl = = J3, jf = .1^ = J4 and = jI = Jl, j1 = = j2, = = J3, j| - J2" = J 4- This 
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implies: 



oiii«n3]ibi) - (jiiianiiibi) =no3,j3)5j3..3+i: (J4ii«n4]iijt) = (J2ii«n2]ib2> =n(j4,j4)%„,,+. 



Using the U(l) identifications and ((83|) . (|84l) . ([87|) . ([88]) and ([90|l and merging all these equations 
carefully, we get: 

{jabcd\T^fUabcd\jabcd) — —S^a VjaS^a ^aS^a (5 b -2b6^b -^bS^c ^c6jc -AC^Aa l^a (Sd -id^ zd ^Jd ( — 1 ^ 

Wuucui uucu ijuut-u/ J3 J3 J4 J4 Jl2'J12 Jl Jl J4 '^4 Jl >Jl J2 J2 J12J12 J2 '■J2 Js J3 ^ ' 

(_i).f+.H.2^+.:?n0u^)n(j2j^)n03j3)n(j4i;)n(ji^jt^ 



6 .-fc MTf-.d \ 1 ■?! 2 1 j J12 J12 2 



H ji & J I Jl jl J4 




■d -d ill nd Ad I 
J12 J12 2 I J -'12 -'12 2 

is js i2 Ji J4 J4 if 



(91) 



In (HH), n(a,6) = (-l)'^+''+5n(a,6). Note that n(a,6) are symmetric n(a,6) = n(6, a) and real. We 
have ignored the 16 S functions Y[i=i ('^Jiji+j coming from the reduced matrix elements 

in (|79|) as they are already contained in the six 6j symbols in ([9T|) . The above d=2 loop dynamics and 
it's generalization to arbitrary dimensions are discussed in detail in sections (3.1) and (3.2). 
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